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Models of the universe with arbitrary (non gravitational) interaction between the components of 
the cosmic fluid: the phantom energy and the background, are investigated. A general form of the 
interaction that is inspired in scalar-tensor theories of gravity is considered. No specific model for 
the phantom fluid is assumed. We concentrate our investigation on solutions that are free of the 
coincidence problem. We found a wide region in the parameter space where the solutions are free 
of the big rip singularity also. Physical arguments, together with arguments based on the analysis 
of the observational evidence, suggest that phantom models without big rip singularity might be 
preferred by Nature. 
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I. INTRODUCTION 

Recently it has been argued that astrophysical ob- 
servations might favor a dark energy (DE) component 
with " supernegative" equation-of-state (EOS) parameter 
uji = Pi/ Pi < —1 Q} Q- where pi is the pressure and pi the 
energy density of the i-th component of the cosmic fluid. 
Sources sharing this property violate the null dominant 
energy condition (NDEC). Otherwise, well-behaved en- 
ergy sources (positive energy density) that violate NDEC 
(necessarily negative pressure), have EOS parameter less 
than minus unity 0. NDEC- violating sources are be- 
ing investigated as possible dark energy (DE) candidates 
and have been called as "phantom" components 0,0 -[27] 
Since NDEC prevents instability of the vacuum or prop- 
agation of energy outside the light cone, then, phan- 
tom models are intrinsically unstable. Nevertheless, if 
thought of as effective field theories (valid up to a given 
momentum cutoff), these models could be phenomeno- 
logically viable J2j. Another very strange property of 
phantom universes is that their entropy is negative |6(. 

To the number of unwanted properties of a phan- 
tom component with "supernegative" EOS parameter 
Uph — Pph/Pph < —1, we add the fact that its en- 
ergy density p p h increases in an expanding universe. [28] 
This property ultimately leads to a catastrophic (future) 
big rip singularity |8( that is characterized by divergences 
in the scale factor a ] the Hubble parameter H and its 
time-derivative £T[l(j- Although other kinds of singu- 
larity might occur in scenarios with phantom energy, in 

the present paper we are interested only in the big rip 

kind of singularity. [29] This singularity is at a finite 

crcncc 

[30] Although experimental tests in the solar system impose severe 

restrictions on the possibility of n on-m inimal coupling between 
* Electronic address: rcurbelofauclv.edu.cuitamcfauclv.edu.cuasraelfauclv.eM.AA 1 ^ and ordinary matter fluids [H, due to the unknown na- 
ture of the DM as part of the background, it is possible to have 
additional (non gravitational) interactions between the dark en- 
ergy component and dark matter, without conflict with the ex- 
perimental data. It should be pointed out, however, that when 
the stability of DE potentials in quintessence models is consid- 
ered, the coupling dark matter-dark energy is troublesome I 111 . 
Perhaps, the argument might be extended to phantom models. 



amount of proper time into the future but, before it is 
reached, the phantom energy will rip apart all bound 
structures, including molecules, atoms and nuclei. To 
avoid this catastrophic event (also called "cosmic dooms- 
day"), some models and/or mechanisms have been in- 
voked. In Ref.0], for instance, it has been shown that 
this singularity in the future of the cosmic evolution 
might be avoided or, at least, made milder if quantum 
effects are taken into consideration. Instead, a suitable 
perturbation of de Sitter equation of state can also lead to 
classical evolution free of the big rip 0] . Gravitational 
back reaction effects an d scalar fields with negative 
term with self interaction potentials with 
have also been considered . 
Another way to avoid the unwanted big rip singularity 
is to allow for a suitable interaction between the phan- 
tom energy and the background (DM) fluid[illlT7j]. [30] 
In effect, if there is transfer of energy from the phantom 
component to the background fluid, it is possible to ar- 
range the free parameters of the model, in such a way 
that the energy densities of both components decrease 
with time, thus avoiding the big rip^(|. Models with in- 
teraction between the phantom and the DM components 
are also appealing since the coincidence problem (why 
the energy densities of dark matter and dark energy are 
of the same order precisely at pres ent?) can be solved or, 
at least, smoothed out [la. Il7l l20| . 

Aim of the present paper is, precisely, to study models 
with interaction between the phantom (DE) and back- 
ground components of the cosmological fluid but, un- 



[27] That NDEC-violating sources can occur has been argued decades 



ago. 



for instance, in reference [fj 



[28] Alternatives to phantom models to account for supernegative 
EOS parameter have been considered also. See, for instance, 
references 

[29] A detailed study of the kinds of singularity might occur in phan- 
tom scenarios (including the big rip) has been the target of ref- 
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like the phenomenological approach followed in other 
cases to specify the interaction term (see references 
[Ifl Il7j). we start with a general form of the interac- 
tion that is inspired in scalar-tensor (ST) theories of 
gravity. We concentrate our investigation on solutions 
that are free of the coincidence problem. In this pa- 
per no specific model for the phantom energy is as- 
sumed. A flat Friedmann-Robertson- Walker (FRW) uni- 
verse (ds 2 = —dt 2 + a(t) 2 Sijdx l dx J ; ; i, j — 1,2,3) is con- 
sidered that is filled with a mixture of two interacting flu- 
ids: the background (mainly DM) with a linear equation 
of state p m = Lo m p m , Lu rn = const., [31] and the phantom 
fluid with EOS parameter uj p h = p p h/pph < ~ !• Addi- 
tionally we consider the " ratio function" : 



equation for each fluid shows a source (interaction) term: 



Pm 

p P h 



between the energy densities of both fluids. It is, in fact, 
no more than a useful parametrization. We assume that 
r can be written as function of the scale factor. 

The paper has been organized in the following way: In 
Section 2 we develop a general formalism that is adequate 
to study the kind of coupling between phantom energy 
and the background fluid (inspired in ST theories as said) 
we want to investigate. Models of references |l6|,LL3 can 
be treated as particular cases. The way one can escape 
one of the most undesirable features of any DE model: 
the coincidence problem, is discussed in Section 3 for the 
case of interest, where the DE is modelled by a phantom 
fluid. In this section it is also briefly discussed under 
which conditions the big rip singularity might be avoided. 
With the help of the formalism developed before, in Sec- 
tion 4, particular models with interaction between the 
components of the mixture are investigated where the 
coincidence problem is solved. The cases with constant 
and dynamical EOS parameter are studied separately. It 
is found that, in all cases, there is a wide range in pa- 
rameter space, where the solutions are free of the big 
rip singularity also. It is argued that observational data 
seems to favor phantom models without big rip singular- 
ity. This conclusion is reinforced by a physical argument 
related to the possibility that the interaction term is al- 
ways bounded. Conclusions are given in Section 5. 



II. GENERAL FORMALISM 

Since there is exchange of energy between the phantom 
and the background fluids, the energy is not conserved 
separately for each component. Instead, the continuity 



[31] Although baryons, which should be left uncoupled or weakly cou- 
pled to be consistent with observations, are a non vanishing but 
small part of the background, these are not being considered 
here for sake of simplicity. Also there are suggestive arguments 
showing that observational bounds on the "fifth " fo rce do not 
necessarily imply decoupling of baryons from DE 21]. 



p m + 3H(p m +p m ) = Q, 



pph + 3H(p ph + p ph ) = -Q, 



(2) 



(3) 



where the dot accounts for derivative with respect to the 
cosmic time and Q is the interaction term. Note that 
the total energy density p T = p m + p ph (p T = p m + P P h) 
is indeed conserved: pr + 3H(pt + Pt) — 0. To specify 
the general form of the interaction term, we look at a 
scalar-tensor theory of gravity where the matter degrees 
of freedom and the scalar field are coupled in the action 



(1) through the scalar-tensor metric x{4>Y 



9ab _ 



s ST = J d 4 xV\F\{§ - \m) 2 



(4) 



where x(0) -2 is the coupling function, C m is the matter 
Lagrangian and p is the collective name for the matter 
degrees of freedom. It can be shown that, in terms of 
the coupling function x(</0> the interaction term Q in 
equations @ and J3J), can be written in the following 
form: 



Q = p m H[a 



da 



(5) 



where we have introduced the following "reduced" no- 
tation: x(a) = x(a)^ 3Um ~ 1 ^ 2 and it has been assumed 
that the coupling can be written as a function of the 
scale factor x = x( a )-[32] This is the general form of the 
interaction we consider in the present paper. [33] Com- 
paring this with other interaction terms in the bibliog- 
raphy, one can obtain the functional form of the cou- 
pling function x m each case. In Ref.^iJ, for instance, 
Q = 3Hc 2 (p ph + p rn ) = 3c 2 Hp m (r + l)/r, where c 2 de- 
notes the transfer strength. If one compares this expres- 
sion with one obtains the following coupling function: 



X(o) = Xo e 



(6) 



where xo is an arbitrary integration constant. If 
c 2 = Cq = const, and r = ro = const., then x = 



[32] If Sst represents Brans-Dicke theory formulated in 
the Einstein frame, then the coupling function x{<t>) = 
Xo ex P (~4>/ + 3/2). The dynamics of such a theory but, for 
a standard scalar field with exponential self-interaction potential 
has been studied, for instan ce, in Ref . |23j . For a phantom scalar 
field it has been studied in |24| . 

[33] We recall that we do not consider any specific model for the 
phantom field so, the scalar-tensor theory given by the action 
Sst serves just as inspiration for specifying the general form of 
the interaction term Q. 
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Xo a 3c °( ro+1 '/ r °. Another example is the interaction term 
in Ref. 0]: Q = 5Hp mi where 5 is a dimensionless cou- 
pling function. It is related with the coupling function x 
(Eqn. J5J) through: 



X(a) = Xo ef "" s , 



(7) 



and, for S — Sq; x(a) — Xo a<5 °- 

If one substitutes (J3J) in J2Jl, then the last equation can 
be integrated: 



Pm — Pm,0 a Xi 



(8) 



where p m ,o is an arbitrary integration constant. 

If one considers Eqn. Q then, Eqn.JJJ can be ar- 
ranged in the form: p p h/p P h + 3(w p /, — l)H — (3tu m — 
l)rH[ad(\nx~ 1 ^ 2 )/da] or, after integrating: 



Pph = P P /i,oexp{- 



da , (iflnv)^ , „ 

— l + 3uv,+ra-V^}, (9) 
a da 



where p p h.o is another integration constant. Using equa- 
tions (J8J, Q and @, it can be obtained an equation 
relating the coupling function x, the phantom EOS pa- 
rameter oj p h and the "ratio function" r: 



X(«) = Xo(— — r) 
r + 1 



(10) 



where, as before, xo is an integration constant. There- 
fore, by the knowledge of w p j, = uj p h(a) and r = r(a), 
one can describe the dynamics of the model under study. 
Actually, if u p h and r are given as functions of the scale 
factor, then one can integrate in Eqn. (|10|l to obtain 
X — x(fl) and, consequently, p m = p m (a) is given through 
Eqn. (JSJ). The energy density of the phantom field can 
be computed through either relationship or Be- 
sides, the Friedmann equation can be rewritten in terms 
of r and one of the energy densities: 



3H 2 = p ph (l + r) = p m ( ), 



(11) 



so, the Hubble parameter H = H(a) is also known. Due 
to Eqn. , the (dimensionless) density parameters f2j = 
Pi/3H 2 (£l m + Q P h = 1) are given in terms of only r: 



1 r 



l + r' 



1 + r' 



(12) 



respectively. It is useful to rewrite Eon. l|10|l . alterna- 
tively, in terms of Q p h and Lo p h- 



X(a) = Xo (1 - n ph ) e-^ *^-"») n ". 



(13) 



Another useful parameter (to judge whether or not the 
expansion if accelerated), is the deceleration parameter 
q = —(1+H/H 2 ), that is given by the following equation: 



q = -i 



3 rWph 



1 



l)r 



l + r 



(14) 



or, in terms of Q p h and to p h' 

q = - [1 + 3uj m + 3(u> ph - u) m )Vl ph ] . (15) 

When the EOS parameter of the phantom field is a 
constant ui p h = —(1 + £ 2 ) < — 1, £ 2 £ R + , we need to 
specify only the behavior of the ratio function r — r(a) 
or, equivalently, of the phantom density parameter £l p h = 
^ P h(ci). In this case the equation for determining of the 
coupling function x <|ll)H or, alternatively, (|13|l . can be 
written in either form: 

x(a)=Xo(^-)e 3 ^ 2+ ^f^, (16) 
r + 1 



X (a)=X0 (l-n ph )e 3 ^ 2 +^I^. 



(17) 



The deceleration parameter is given, in this case, by q = 
-1 + (3/2) [(r - t 2 )/(r + 1)] so, to obtain accelerated 
expansion r < 2 + 3£ 2 . 



III. AVOIDING THE COINCIDENCE PROBLEM 

In this section we discuss the way in which one of the 
most undesirable features of any DE model (including the 
phantom) can be avoided or evaded. For completeness 
we briefly discuss also the conditions under which the 
doomsday event can be evaded. 



A. How to avoid the coincidence problem? 

It is interesting to discuss, in the general case, un- 
der which conditions the coincidence problem might be 
avoided in models with interaction among the compo- 
nents in the mixture. In this sense one expects that a 
regime with simultaneous non zero values of the den- 
sity parameters of the interacting components is a criti- 
cal point of the corresponding dynamical system, so the 
system lives in this state for a sufficiently long period of 
time and, hence, the coincidence does not arises. 

In the present case where we have a two-component 
fluid: DM plus DE, one should look for models where 
the dimensionless energy density of the background (DM) 
and of the DE are simultaneously non vanishing dur- 
ing a (cosmologically) long period of time. Otherwise, 
we are interested in ratio functions r(a) that approach 
a constant value tq ~ 1 at late times. A precedent of 
this idea but for a phenomenologically chosen interac- 
tion term can be found, for instance, in Ref. [T^ . where 
a scalar field model of phantom energy has been ex- 
plored (a self-interacting scalar field with negative ki- 
netic energy term p p h — — 2 /2 + V{4>)). In this case, 
as already noted, the interaction term is of the form 
Q = 3Hc 2 (p ph + Pm ) = 3c 2 Hp m {r + l)/r. For this 
kind of interaction, stability of models with constant ra- 
tio tq = p m / p p h (and constant EOS parameter) yields to 
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very interesting results. In effect, it has been shown in 
[l6| that, if r < 1 (phantom-dominated scaling solution) , 
the model is stable. The solution is unstable whenever 
ro > 1 (matter-dominated scaling solution). One can 
trace the evolution of the model that evolves from the 
unstable (matter-dominated) scaling solution to the sta- 
ble (phantom-dominated) one. As a consequence, this 
kind of interaction indicates a phenomenological solution 
to the coincidence problem Actually, once the uni- 
verse reaches the stable phantom-dominated state with 
constant ratio p m /Pph = ^m/^ph — r o> h will live in 
this state for a very long time. Then, it is not a coinci- 
dence to live in this long-living state where, since ro ~ 1, 

Pm ~ Pph • 

The above stability study can be extended to situa- 
tions where no specific model for the phantom energy is 
assumed (like in the case of interest in the present paper) 
and where a general interaction term Q is considered (in 
our case Q is given in Eqn.JSJ)), if one introduces the 
following dimensionless auantities|ll|: 



PK 

3iJ 2 '' 



V 



Pv 
3H 2 



(18) 




where the "kinetic" px = (pph + Pph)/2 and "potential" 
pv = [pph —Pph)/% terms have been considered. Besides, 
if one where to consider FRW universe with spatial cur- 
vature fc, an additional variable: 



(19) 



a 2 H 2 



is needed. The governing equations can be written as the 
three-dimensional autonomous system: 



FIG. 1: The phase portraits (x,y) = (^, f^) for the 
flat space model with constant ratio ro and EOS parameter 
Uph.o = — (1 + £ 2 ) are shown for background dust (oj m = 0). 
In the upper part of the figure we have considered ro < £ 2 
(ro = 1, £ 2 = 3). The point (—3/2,5/2) is an attractor. In 
the lower part of the figure ro > £ 2 (ro = 1, £ 2 = 1/3). 
All trajectories in phase space diverge from an unstable node 
(third critical point (x,y) = (-1/12, 7/12)). The first critical 
point is a saddle. In any case the model has a late time phan- 
tom dominated attractor solution (the second critical point 
(x,y) = (-1/6, 7/6)). 



-(l+3^)(3 
dpph 



Q 



2H PK 



)x 



+3x[2x + -z + (w m + 1)(1 - x - y — z)], 
_(l_^) [3 z + n 



dpph ' 



+3y[2x +-^ z + ( u rn + - x 
2 

-2z + 3z[2x + -z + (oj m + 1) 

O 



(1 



v- z )] 



(20) 



together with the constrain Q m = 1 — x — y — z, where 
the comma denotes derivative with respect to the variable 
N = In a. The parameter range is restricted to be < 
x + y + z < 1. In the remaining part of this section, for 
simplicity, we assume that to m = 0, i. e., the background 
fluid is dust. [34] 



[34] The details of the stability study will be given in a separate 
publication, where other kinds of interaction are considered also. 



1. Flat FRW (k = 0j 

Let us consider first the spatially flat case k = => 
z = 0. Consider, besides, a constant ratio r = ro and a 
constant EOS parameter u! p h = uiph.o = —(1 + £ 2 ) < — 1. 
According to Ijlbl , one obtains a coupling function of the 

form: £(a) = xo (r /(r + 1)) a^ 1+ ^^ ro+1 \ If one 
substitutes this expression for the coupling function into 
Eqn. ©, then: 



Q = — — r~^ H Pn 



ro 



1 



(21) 



The fixed points of the system H2l)fl can be found once 
the interaction function Q given in Ean. H21|) and the ex- 
pression dpph/dpph = ojphfi = —(1 + £ 2 ) are substituted 
into (J20J- The fixed points of the system (x', y') — (0, 0) 
are: 1) (x,y) = (0,1), 2) (x,y) - (-£ 2 /2, 1 + £ 2 /2) and 
3) (x,y) = (-£ 2 /[2(r + l)],(2+a/[2(r + l)]). Of these, 
the first point is a saddle if r > £ 2 and an unstable node 
if r < £ 2 . The second critical point is always a stable 
node. In this case fl p h = x + y = 1 so, there is phantom 
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(dark energy) dominance. The third point is an unstable 
node whenever r > £ 2 and a saddle otherwise. In this 
case ftph = l/(ro+l) and Q m — ro/(ro+l) are simultane- 
ously non vanishing so, it corresponds to matter-scaling 
solution. In Fig. ^ the convergence of different initial 
conditions to the attractor (phantom dominated) solu- 
tion is shown in the phase space (x, y) (flat space case) 
for given values of the parameters ro and £ 2 . If ro > £ 2 
all trajectories in phase space diverge from an unstable 
node (third critical point). The first critical point is a 
saddle. In any case the model has a late-time phantom 
dominated attractor solution (the second critical point). 



2. FRW with spatial curvature (k ^ 0) 

The critical points (x' = 0,y' = 0, z' = 0) are: 
1) (x,y,z) = (0,1,0) (de Sitter space), 2) (x,y,z) = 
(-£ 2 /2, 1 + e/2, 0), 3) (x, y, z) = (-£ 2 /[2(r + 1)], (2 + 
£ 2 )/[2(r + l)],0) and 4) (x,y,z) = (0,0,1) (Milne space). 
The first point is always a saddle, meanwhile, the second 
critical point is always a stable node (a sink). The third 
critical point is a saddle whenever ro < £ 2 and an un- 
stable node (a source) otherwise. This point corresponds 
to matter-scaling solution where both Q m and tt p h are 
simultaneously non negligible. The fourth critical point 
(Milne space) corresponds to curved (k ^ 0) vacuum 
space (f2 m = £l p h — 0). This point is an unstable node if 
ro < 2 + 3£ 2 and a saddle whenever ro > 2 + 3£ 2 . Note 
that, even if the evolution of the universe begins with 
initial conditions close to the fourth critical pint (curved 
flat space), the model drives the subsequent evolution of 
the universe to the stable node corresponding to a flat 
universe dominated by the phantom component. 

In summary, the way to evade the coincidence prob- 
lem when an interaction term of the general form @ is 
considered, is then clear: One should look for models 
where, whenever expansion proceeds, both u p h and the 
ratio r = p m /Pph — (l/^ph) — 1, approach constant val- 
ues uiph.o and ro ~ 1, respectively. Then, if trajectories in 
phase space move close to the third critical point (matter- 
scaling solution), these trajectories will stay close to this 
point for an extended (arbitrarily long) period of time, 
before subsequently evolving into the phantom domi- 
nated solution (second critical point) that is a late-time 
attractor of the model. [35] 



B. How to avoid the big rip? 

If we assume expansion (the scale factor grows up) 
then, to avoid the future big rip singularity, it is nec- 
essary that the background fluid energy density p m and 
the phantom pressure p p h and energy density p p h, were 



bounded into the future. Once the "ratio function" 
r (Eqn.Q) and the parametric EOS for the phantom 
fluid (p p h — u P hP P h) are considered, the above require- 
ment is translated into the requirement that p p h, 0J p hP P h 
and rp p h were bounded into the future. Let us con- 
sider, for illustration, the case with constant EOS: to p h — 
Uph.o = — (l + £ 2 ) < — 1, and the following ratio function: 



r+ E3- 



Integrating in the exponent 



in Ean. H16|) one obtains the following expression for the 
coupling function: 

x(a) = Xor (a k + r Q )- 1+3 ( 1+ ^/ k . (22) 

The expression for the DM and phantom energy densities 
are then given through equations (see equations (JSJ) and 



and 



Pm=Pm,a a- 3 (a fc +ro)- 1+3 ( 1+ « 2 )/ fc , 



Pph = P P h,0 « 



-3+k 



(a" 



ro) 



-l+3(l-K 2 )//c 



(23) 



(24) 



where p mfi = p m ,oXo r Q and p phfi = PmfiXo- We see 
that, for a ^ rj' k , p m oc a~ fc + 3 ? i while p p h oc a 3 ^ . 
For £ 2 < fc/3 the DM energy density p m decreases with 
the expansion of the universe. However, the phantom 
energy density always increases with the expansion and, 
consequently, the big rip singularity is unavoidable in this 
example. The conclusion is that, for a phantom fluid 
with a constant EOS uo p h — u p h,o and a power-law ratio 
of phantom to background (dust) fluid energy densities 
Pm/p P h oc a~ k , the big rip singularity is unavoidable. 

For a constant ratio r = ro and a constant EOS (u> p h — 
u p h.o = —(1 + £ 2 ) < —1) UMi direct integration in Eqn. 
(11611 yields to the following coupling function: 



X 



ro 



( a ) = Xo ( — r) a -o+i 

ro + 1 



(25) 



so, the background and phantom energy densities in Eqn. 
||SJ arc given by: 



P P h — 7*0 P n 



Pm,0 -3(2^£-) 

— '— a K r o+ 1 ', 
ro 



(26) 



where p m ,o = Pm,oXo [ r o/(^o + l)]- We see that, whenever 
there is expansion of the universe, for ro > £ 2 , there is 
no big rip in the future of the cosmic evolution, since 
p ph oc r p ph oc uJ p h,oPph oc a~ 3(r °^")/( ro+1 ) are bounded 
into the future. In order to get accelerated expansion, 
one should have ro < 2 + 3£ 2 so, the ratio of background 
(dust) energy density to phantom energy density should 
be in the range £ 2 < r < 2 + 31; 2 . 



IV. MODELS THAT ARE FREE OF THE 
COINCIDENCE PROBLEM 



[35] This argument is due to Alan Coley (private communication). 



In the present section we look for a model with an 
appropriated coupling function \ that makes possible to 



6 
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redshif t 



FIG. 2: The plot of energy density parameters Q m (dashed 
curve) and Q. p h (thin solid curve) as well as of the ratio func- 
tion r (thick solid curve) vs redshift z, is shown for the model 
with constant EOS parameter. For simplicity background 
dust (w m — 0) was considered. The values of the free pa- 
rameters are: m = 12, C = 0.02 and tv p h,a ~ —1.1. The 
ratio function r approaches a constant value at negative red- 
shift ro = 0.23. At present (z = 0) r = 0.25, which means 
that we are already in the long living matter-scaling state. 
The background energy density dominates the early stages of 
the cosmic evolution. At z ~ 0.33 both density parameters 
equate and, since then, the phantom component dominates 
the dynamics of the expansion. A very similar behavior is 
obtained for fi m , Q p h and r in the case with dynamical EOS 
parameter (Case B). 



avoid the coincidence problem. We study separately the 
cases with a constant EOS parameter u> p h — w P h,o — 



-(1- 
ctcr. 



-£ 2 ) < — 1; £ 2 6 R + and with variable EOS param- 



A. Constant EOS parameter uj p h = u) p h,o 

In order to have models where both f2 m and £l p h are 
simultaneously non-vanishing during a long interval of 
cosmological time (no coincidence issue) and that, at the 
same time, are consistent with observational evidence on 
a matter dominated period with decelerated expansion, 
lasting until recently (redshift z ~ 0.39 2]), we choose the 
following dimensionless density parameter for the phan- 
tom component: 



%h = ~s ( 



B a ri 



C 



(27) 



where m, B, and C are arbitrary constant parameters. 
This kind of density parameter function can be arranged 
to fit the relevant observational data, by properly choos- 
ing the free parameters (m, B, C) (m controls the curva- 
ture of the curve £l p h(z), meanwhile C controls the point 
at which £l p h(z eq ) = ^ m {z eq )). In this sense, the model 
under study could serve as an adequate, observationally 
testable, model of the universe. The choice of Q p h in 



equation 12 7|) . when substituted in Eqn. 1|17[) yields to 

(1 - § )a m + C 



xip) = Xo 



(a m + C) 1 - 3 ( 



U) m —UJ p h 



o)/B- 



(28) 



In consequence, equations © and lead to the follow- 
ing expressions determining the background energy den- 
sity and the energy density of the phantom component 
respectively; 



Pm — Pm,o a 



-3 



(1 - %)a m + C 



(a ,n + C , ) 1 -3(w m -w phi0 )/-B ' 
where p mfi = p m . Xo/B, p = p m ,o£o and 



(29) 



P P h = P P h,o a m - 3 (a m + C )-i+3(<*»-u>,>.o)/b (30) 

where p phfi = m p m , Xo/B. 

It is worth noting that, if uj m — B/m < Lo p h.a{< — 1), 
there is no big rip singularity in the future of the cosmic 
evolution as described by the present model. In effect, 
for large a 3> 1, the following magnitudes 



p P h oc r p ph oc uj phfi p ph 
•y „-3(l-m(w m -w p h,o)/B) 



(31) 



are bounded into the future. In particular, the phan- 
tom energy density is a decreasing function of the scale 
factor. Note that the ratio function r = Q m /£l p h ap- 
proaches the constant value (B/m) — 1 so both fi m and 
fl p h are simultaneously non vanishing. As explained in 
the former section, this means that the universe is in the 
third critical point (matter-scaling solution) which is a 
saddle for r < £ 2 . In consequence, the universe will 
evolve close to this critical point for a sufficiently long 
period of time and then, the model approaches a stable 
dark energy dominated regime (second critical point). 

In Fig. 0we show the behavior of the density parame- 
ters Q m and Q p h and of the ratio function r as functions 
of the redshift (for simplicity we have considered back- 
ground dust: uj m — 0). The following values of the free 
parameters: m = 12, u> p h,o = — 1.1, C = 0.02 have 
been chosen, so that it is a big rip free solution. Be- 
sides we set the DE density parameter at present epoch 
Qph(z = 0) = 0.8 so, the following relationship be- 
tween the free constant parameters B, m and C should 
be valid: B = 1.25[m/(l + C)]. In consequence, for 
the values of m and C given above: B/m = 1.22 so, 
ujph.o = —1.1 > —1.22, meaning that there is no big 
rip singularity in the future of the cosmic evolution in 
this model (see the condition —B/m < uj p h,o under Eqn. 
(|30|l ). In Fig. the plot of the deceleration param- 
eter q vs redshift is shown for three different values of 
the EOS parameter u p h,o = —1-1, —1.5 and —3 respec- 
tively. It is a nice result that the solution without big rip 
(thick solid curve) is preferred by the observational evi- 
dence, since, according to a model-independent analysis 
of SNIa data[25| , the mean value of the present value of 
the deceleration parameter is < qo >« —0.76. 
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FIG. 3: The deceleration parameter g is plotted as function 
of the redshift z for the Model A (constant EOS parame- 
ter), for three different values of = —1.1, —1.5 and 
—3 (thick solid, thin solid and dashed curves respectively). 
Note that the transition from decelerated (positive g) into 
accelerated (negative g) expansion occurs at z ~ 0.4 for the 
value Uphfi = — 1.1. It is apparent that this is, precisely, the 
observationally favored solution, since, according to a model- 
independent analysis of SNIa data, q(z = 0) ~ — 0.76|25|). In 
consequence, solutions without big rip are preferred by obser- 
vations. 



The way the coincidence problem is avoided in this 
model is clear from FigEl also. The ratio between the 
energy density parameters of the background and of the 
DE approaches a constant value r = 0.23 for negative 
redshift so, both fi m and Q, p h are non negligible. As 
already said, this is a critical point (third critical point 
in the stability study of section 3) and the universe can 
live for a long time in this state. Moreover, since at 
present (z = 0); r(z = 0) = 0.25, we can conclude that 
the universe is already in this long living state. 



B. Dynamical EOS parameter 

We now consider a dynamical EOS parameter u> p h = 
u> v h(a). In this case we should give as input the functions 
Q p h(a) and u> p h{a). In consequence, one could integrate 
in the exponent in Ean. l|13l) so that the dynamics of the 
model is completely specified. In order to assure avoid- 
ing of the coincidence problem, let us consider the same 
phantom energy density parameter as in the former sub- 
section (Eqn. (|27|l ). but rewritten in a simpler form: 



a a" 



(32) 



where m, a and f3 are non negative (constant) free param- 
eters. To choose the function uj p h(a) one should take into 
account the following facts: i) at high redshift and until 
recently [z ~ 0.39 ± 0.03 0) the expansion was deceler- 
ated (positive deceleration parameter q) and, since then, 
the expansion accelerates (negative q), ii) for negative 
redshifts the EOS parameter approaches a constant value 
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FIG. 4: The plot of the energy densities of DM p m (upper part 
of the figure), phantom component p p h (middle part) and of 
the Hubble parameter H vs redshift, is shown for three values 
of the constant parameter ui p h,o- —1.1 (thick solid curve), 
— 1.5 (thin solid curve) and —3 (dashed curve) for the model 
with dynamical EOS parameter. It is apparent that, only in 
the first case the model is big rip free. Note that, due to the 
values chosen for the free parameters, in particular m = 12, 
the phantom energy density vanishes as one approaches the 
initial big bang singularity. This is true whenever m > 3 (see 
Eqn. JUJ). 



Wph,Q more negative than —1; uj p h,o — — (l+£ 2 ), £ 2 £ R + ■ 
We consider, additionally, the product (uj p h — u m ) Q p h to 
be a not much complex function, so that the integral in 
the exponent in Eqn. 11311 could be taken analytically. A 
function that fulfills all of the above mentioned require- 
ments is the following: 

( \ j. (a m + 0Ka m -5) 

u p h(a) — Lo m + LUphQ , (66) 

a zm + o 

where 5 is another free parameter. The parameter m 
controls the curvature of the density parameter function 
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FIG. 5: The interaction term Q is plotted vs redshift z. 
Note that for 0.2 ;$ z ;S 1.7, Q is negative, meaning that 
there is transfer of energy from the DM into DE component. 
For higher z the model proceeds without interaction. Q is 
bounded into the future only for the big rip free case (thick 
solid curve). 



Qph(z), while d controls the point of equality £l p h{Zeq) 
VL m (z eq ). The resulting coupling function is 



X(a) = Xo (1 - tt ph ) (a 2m + 5)-fe"**.o 

; ^2: arctan(^) 



x e 



(34) 



Therefore, the background and phantom anergy densities 
are given by the following expressions: 



p m (a) = p m ,o m ' R a 3 ^ +1 > 

a m + p 

,2m , X \—iSrU>„ h . a „^^arctan(^-) 



(35) 



PphW = Pph,0 



m— 3(uJ m + l) 



(a 2m +6)-^^° 

x e "> y ^ ' , (36) 



where the constant p p h,a = en p m ,o- At late times 
(large a), the functions r p p h <x u> p h p P h oa p p h <x 
a -3(u> m +i+au phi0 ) _ This meang that! w henever -(1 + 

w OT )/a < w p / l) o(< —1), these functions are bounded into 
the future and, consequently, there is no big rip in the 
future of the cosmic evolution in the model under study. 

In FigQ] the plot of the energy densities of DM p m 
(upper part of the figure) , phantom component p p h (mid- 
dle part) and of the Hubble parameter H vs redshift, is 
shown for three values of the constant parameter u> p h,o- 
—1.1 (thick solid curve), —1.5 (thin solid curve) and —3 
(dashed curve). The following values of the free param- 
eters: m = 12, j3 = 0.03, and 5 = 3 10 -4 have been 
chosen. We assume Q m (z = 0) = 0.3j2|J so, the fol- 
lowing relationship between a and /3 should take place: 
a = 0.7(/3 + 1). It is apparent that, only in the first 
case (<jJ p h,o = ~ 1-2) the model is big rip free. Actu- 
ally, as already noted, the necessary requirement for ab- 
sence of big rip is, for the chosen set of parameters, 
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FIG. 6: In the upper and middle part of the figure the EOS 
parameter uj p h and the deceleration parameter q are plotted 
as functions of the redshift z (Model B), for three values of 
the constant parameter uj p h,o = — 1.1, —1.5 and —3 (thick 
solid, thin solid and dashed curves respectively). It is worth 
noting that at high redshift uj p h approaches a constant value 
u p h =0. It is a minimum at z = 0.35. The dark energy 
became a phantom fluid only recently, at z ~ 0.45. Note that 
the transition from decelerated (positive q) into accelerated 
(negative q) expansion occurs at z ~ 0.4 for the value u> p h,o = 
— 1.1. This last (big rip free) solution, seems to be preferred by 
observations since, according to Ref. [2^, q(z = 0) » —0.76. 
In the lower part of the figure the ratio function is plotted 
for three values of the parameter m: 12 (thick solid curve), 
8 (thin solid curve) and 4 (dashed curve) respectively. In 
any case, at present, the universe is already in the state with 
constant ratio r. 



LUphfl > —1/a = —1.387. It is noticeable that, since 
we chose m = 12 > 3, the DE energy density dimin- 
ishes as one goes back into the past (it is a maximum at 
present) . In all cases the interaction term is negative dur- 
ing a given period of time in the past (see Fig[SJl, meaning 
that the DM transferred energy to the DE component. 
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At higher redshift the evolution proceeds without inter- 
action. Worth noting that only in the singularity free 
case (uiph,o = —1-2) the interaction term is bounded into 
the future. 

A physical argument against doomsday event can be 
based on the following analysis. At late time, the inter- 
action term Q ~ —3a u> p h,a p m H. This means that, 
only in the case when the evolution proceeds without big 
rip (pm and H are bounded), the interaction between 
the components of the cosmic fluid does not become un- 
physically large at finite time into the future. Besides, 
for the big rip free the expansion proceeds into 

the future, the Hubble parameter approaches a negligible 
value, meaning that the ratio Q/p m is small. 

In the model the transition from decelerated into ac- 
celerated expansion takes place at z ~ 0.4 (see Fig. EJ. 
At late times the deceleration parameter approaches the 
constant value qo f» —0.8. Therefore, as in the case with 
constant EOS parameter (Model A), the observational 
evidence (as argued in Ref. [25|. q(z = 0) « —0.76) seems 
to favor solutions without big rip singularity. It is inter- 
esting to note that, in this model, the dark energy compo- 
nent behaved like dust (ujph — 0) at early times (high z), 
i. e., it behaved like an "ordinary" fluid with attractive 
gravity and, just until very recently (z ~ 0.45), it has not 
became a phantom fluid. At present, the universe is al- 
ready in a state characterized by constant uj p h,o ~ —1.2 
and r = 0.48, meaning that fl m ,o = r /(r + 1) and 
&ph,o = l/( r o + 1) ar< 3 simultaneously non vanishing. As 
shown before, this is a critical point of the corresponding 
dynamical system so, probably, the universe will stay in 
this state for a long time and, consequently, the coinci- 
dence problem does not arise. 



V. CONCLUSIONS 

In this paper we have investigated models with addi- 
tional (non gravitational) interaction between the phan- 
tom component and the background. This kind of in- 
teraction is justified if the interacting components are 
of unknown nature, as it is the case for the DM and 
the DE, the dominant components in the cosmic fluid. 
For the sake of simplicity, baryons, a non vanishing, but 
small component of the background, have not been con- 
sidered in the present study, however, in a more complete 
study these should be included also. A comment should 
be made in this regard: There are suggestive arguments 
showing that observational bounds on the "fifth" force do 
not necessarily imply decoupling of baryons from DE 21]. 

Unlike the phenomenological approach followed in 
other cases to specify the interaction term (see refer- 
ences |Ty, HjJ), we started with a general form of the 
interaction that is inspired in ST theories of gravity (see 
Eqn.||SJ). We have considered no specific model for the 
phantom fluid. Within this context, cosmological evolu- 
tions without coincidence problem have been the target 
of this investigation. We have studied two different mod- 



els: Model A, where the DE EOS parameter is always a 
constant (phantom DE) and Model B, where we consid- 
ered a dynamical DE EOS parameter (the DE became a 
phantom just recently). In these models there is a wide 
range in the parameter space where the cosmic evolution 
is free of the unwanted big rip singularity also. 

The interaction between the different components of 
the cosmic mixture, enables the cosmic evolution to pro- 
ceed without big rip singularity. In effect, if the inter- 
action term is chosen so that there is transfer of energy 
from the phantom component to the background, then it 
is possible to arrange the free parameters of the model 
to account for decreasing energy density of the phantom 
fluid (as well as of the background fluid). A physical ar- 
gument against big rip can be based on the analysis of 
the interaction term Q. In model B, for instance, at late 
time the interaction term Q ~ —3a w p h o Pm H. This 
means that, only whenever p m and H are bounded (no 
big rip), the interaction between the components of the 
cosmic fluid does not become unphysically large at finite 
time into the future. Besides, for the big rip free case, 
as the expansion proceeds into the future, the Hubble 
parameter approaches a negligible value, meaning that 
the ratio Q/p m is small. From the observational per- 
spective it seems that the evolution without doomsday 
event is preferred also. For this purpose we have consid- 
ered model-independent analysis of SNIa data yielding 
to a mean value of the present value of the deceleration 
parameter < qo >« — 0.76|25j. 

Interacting models of dark energy are also useful to 
account for the coincidence problem. In the present 
study, it has been shown that, a solution with simul- 
taneously non vanishing values fl m ,o = ^o/( r o + 1) and 
Qphfi = V( r o + 1) ( r o is the constant ratio of the DM to 
DE densities), is a critical point of the corresponding dy- 
namical system (third critical point in section 3). Hence, 
if trajectories in phase space get close to this point, these 
trajectories will stay there during a sufficiently long pe- 
riod of time. This means that, once the universe is driven 
into the state with simultaneously non vanishing Q m and 
Sl p h, it will live in this state for a (cosmologically) long 
period of time until it is finally attracted into the phan- 
tom dominated phase being a stable node of the models 
(second critical point in section 3). This is the way the 
coincidence problem is solved in the present investiga- 
tion. It seems that we live already in the matter-scaling 
regime (third critical point) and will live in it for long 
time until, finally, the phantom dominates the destiny 
of the cosmic evolution. In any case the evolution of the 
universe will proceed without the risk of any catastrophic 
(doomsday) destiny. 

An obvious limitation of the present study is that only 
two specific models (model A and model B) have been 
considered. However, we think this suffices to illustrate 
how interacting models can be constructed that avoid 
the coincidence problem and , at the same time, fit some 
of the observational evidence. The outcome that big rip 
free models are preferred by some of these observations 
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is just a nice result. 
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